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Abstract 



We report new measurements of the differential and total branching ratios for inclusive B 
decay to D° , D + and D* + and the first measurement of the same quantities for inclusive 
B decay to D*°. Here B is the mixture of Bd and B u from T(4S) decay. Furthermore, 
since more than one charm particle (or antiparticle) of the same kind can be produced in 
B decay, here "inclusive B branching ratio" is used to mean the average number of charm 
particles and their antiparticles of a certain species produced in B decay. We obtain the 
following results (the first error is statistical, the second systematic of this analysis, the 
third is propagated from other measurements): B(B -> D°X) = (0.636 ± 0.014 ± 0.019 ± 
0.018), B(B -> D+X) = (0.235 ± 0.009 ± 0.009 ± 0.024), B(B -> D*°X) = (0.247 ± 0.012 ± 
0.018 ± 0.018), B{B -> D*+X) = (0.239 ± 0.011 ± 0.014 ± 0.009). The following ratio of 
branching ratios is not affected by most of the systematic errors: B(B -> D*°X)/B{B -» 
D* + X) = (1.03±0.07±0.09±0.08). We also report the first measurement of the momentum- 
dependent D*° polarization and a new measurement of the D* + polarization in inclusive 
B decay. Using these measurements and other CLEO results and making some additional 
assumptions, we calculate the average number of c and c quarks produced in B decay to be 
(n c ) = 1.10 ±0.05. 
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I. INTRODUCTION 



Measurements of weak decays of B mesons are essential to testing and understanding the standard model and 
determining the fundamental flavor-mixing parameters. These measurements also provide a unique opportunity to 
examine the short distance behaviour of weak decays JiJ]. Due to their large energy release, long distance corrections 
are expected to be less significant in hadronic B decay than in hadronic decays of charm and strange quarks. However, 
the formation of hadrons in the final state is still poorly understood and is an obstacle to predicting rates and spectra 
for hadronic decays ||[|. 

For several years the small observed value of the B semileptonic decay branching fraction, Bsl(B), seemed to 
disagree with theoretical calculations • Recently, higher order perturbative calculations, taking also into 

account the charm quark mass pO[pj| , come close to reconciling theory with experiment in this respect. However, 
the low value of Bsl(B) implies a larger than naively expected non-leptonic B decay rate ]T^ |. Two mechanisms 
have been proposed and discussed in the literature [P, p^|Jl]] , p^ , p^|Jl4| ] : an enhancement of b — > ccs or of b — > cud. 
The former would give a larger average number of c and c quarks, \n c ), per B decay. Since B -> D° and B -> D + 
transitions (where the D can also be the decay product of charm resonance) account for a large fraction of the charm 
quarks produced in B decay, it is important to measure accurately the branching fractions for these transitions. The 
shape of the momentum spectrum can also be compared to models of hadronic B decay and is sensitive to new 
production mechanisms pUl2|]. In addition, the D and D* inclusive spectra are of interest to future B experiments 
and high energy colliders as Monte Carlo simulations must be constrained to agree with the observed production at 
the T(45). 

In this paper, we describe high statistics measurements of the differential and total branching ratios for the inclusive 
B decay to D and D* mesons, including the first measurement of B(B — ► D*°X). For previous measurements, see 
references |l5| , [l6| . 

We have also measured the D* + and D*° polarizations as a function of the D* momentum. Since there is no 
complete reconstruction of the B final state, we do not distinguish between B° and B + . We produce BB states from 
T(4S) decays, and our generic B is about an even admixture of Bd and B u mesons and their antiparticles . In 
the following, reference to the charge conjugate states is implicit unless explicitly stated. 

It is possible for a B meson to decay to a final state containing two D mesons. By "inclusive B decay branching 
ratios" to a given D species measured in this analysis we mean the "average number of D and D" per B decay. 

After a brief description of the CLEO II apparatus and data sample used, we describe our analysis procedure in 
Sect. III. In Sect. IV- VI we give the B branching ratios and momentum spectra results and analysis details specific to 
the four individual channels. In Sect. VII, we show the results on D* polarization and in Sect. VIII we discuss and 
summarize our results. 



II. DETECTOR AND DATA SAMPLE 



The data sample used in this inclusive B decay analysis was taken at CESR, the Cornell Electron Storage Ring, 
during 1990 - 1994. The sample consists of about 2,020 pb" 1 of integrated luminosity of e + e~ annihilation data 
taken near the peak of the T(4S) resonance and about 959 pb^ 1 just below the open bottom threshold (referred to 
in this paper as "continuum" data). The data used correspond to about 2.166 million BB events. We estimate the 
actual number of BB events by subtracting the number of events in the continuum data after scaling by the ratio of 
the luminosities and correcting for the center-of-mass energy dependence of the continuum annihilation cross-section. 
We also perform minor corrections due to minor differences in CESR running conditions at the two energies. This is 
equivalent to multiplying the integrated luminosity collected near the peak of the T(4S) resonance by the peak cross 
section a(e + e~ — » T(4S)) = (1.072 ± 0.019) nb extracted from our own data. For the purpose of this analysis we 
assume that the T(4S) decays exclusively to BB. 

In the CLEO II detector Jll| charged particles are tracked in a 1.5 Tesla magnetic field through three nested 
coaxial cylindrical drift chambers covering 94% of the solid angle. The innermost chamber is a six-layer straw-tube 
vertex detector of inner radius 4.5 cm with 50 pan position accuracy in the r — (j> plane. It is followed by a ten-layer 
pressurized inner drift chamber with a position accuracy of 100 /im in r — 4>. The main cylindrical drift chamber fllEI] 
contains 51 anode layers, 11 of which are strung at angles to the z-axis progressing from about ±4° to about 7°. It has 
a position accuracy of 120 p.m in r — (f>, it gives a transverse momentum resolution of (Spt/pt) 2 = (0.0015pt) 2 + (0.005) 2 
(with p t in GeV) and a dE/dx resolution (measured on Bhabha scattering events) of 6.5% for particle identification, 
giving good n / K separation up to 700 MeV/c. The outer radius of the main drift chamber is 1 m. Cathode layers 
are located at the inner and outer radii of the ten-layer inner drift chamber, and at the inner and outer radii of the 
main drift chamber, to improve information about the z coordinates (along the beam) of the tracks. Time of flight 
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counters with 154 ps resolution are located outside the drift chambers and provide additional information for particle 
identification (not used in this analysis). 

Photon and ir° detection as well as electron identification use the Csl electromagnetic shower detector pc| . It 
consists of 7800 CsI(Tl) crystals between the time of flight counters and the super-conducting magnet coil in the 
barrel region and, in endcaps, between the drift chamber plates and the magnet pole pieces, altogether covering 95% 
of the solid angle. 

The material in the drift chamber endplates, electronics and cables degrades the performance of the calorimeter in 
the endcaps, especially at the two ends of the barrel region. Photon candidates are restricted to lie in the region of 
the calorimeter covering the angular region | cos#| < 0.707. The energy calibration makes use of Bhabha scattering 
and radiative Bhabhas as well as e + e~ — > 77 reactions and tt° —* 77 decays EL For low multiplicity final states, 
the energy resolution in the barrel portion of the calorimeter is given 

by ^(%) = 0.35/ E - 75 + 1.9- 0.1E, where E 

is the photon energy in GeV. 

Muons can be identified by their penetration in the three 36 cm thick slabs of iron that surround the super- 
conducting coil in an octagonal geometry, and in the iron pole pieces of the magnet 

A. Monte Carlo simulation 

To estimate detection efficiencies we generated Monte Carlo events using the Jetset 7.3 package for continuum 
annihilation events and the CLEO model for T(4S') — > BB decays [M. Separate Monte Carlo datasets were generated 
for the analyses of B — > D°X and of B — > D + X decays. In the Br case, the D + is allowed to decay only into the 
final state we use, K~-k + tv + , while the D~ is allowed to decay according to measured branching fractions of its 
decay modes. About 110,000 BB events and about 35,000 continuum events that contain a D + were generated. An 
analogous procedure is followed in generating D° — * K~tv + , except that the D° is also allowed to decay into K + K~ 
and n + TT~ final states, according to their measured decay branching fractions relative to the K~tt + mode. About 
170,000 BB events and about 150,000 continuum events that contain a D° were generated. The events are then 
processed through a GEANT-based J2f| simulation of the CLEO II detector and reconstructed and analyzed as real 
data. We call these Monte Carlo datasets the "dedicated" Monte Carlo. 

We also used a set of Monte Carlo events produced in a similiar manner, but with all D mesons decaying according 
to a model which incorporates the current knowledge of their decay modes. This statistically independent "generic" 
Monte Carlo (in which the "true" value of the quantities we aim to measure is known a priori) was used to check the 
analysis procedure. This type of consistency check, of course, does not exclude the possibility of a systematic flaw 
that is common to both the "dedicated" and "generic" Monte Carlo data samples. 

The "generic" Monte Carlo was also used for a different purpose. The reconstructed tracks were associated with 
the simulated particles that generated them. We could thus generate high statistics, background free distributions of 
the signals and distributions of various backgrounds. We called these distributions Monte Carlo Tagged (MC-tag for 
short) distributions. 

III. ANALYSIS PROCEDURE 

We reconstruct the D and D* mesons using the following exclusive decay modes that have the largest signal to 
background ratios and for which the Monte Carlo can reliably estimate the reconstruction efficiency: 



D a ^K--K + (1) 

D + -» K~tt + tt + (2) 

D *° _ D \o _» (k-tt + )tv° (3) 

D*+ -► D\ + -» {K-tx + )tx + (4) 

D* + -> D + ir° -» {K--K + TT + )ir° (5) 



In this section, we describe the selection criteria and general procedures used in the analysis. In the sections which 
follow the details of the procedures for decay modes of the D°, D + , D*° and D* + will be given. 

A. Selection Criteria 

The continuum background is suppressed by excluding candidates from jet-like events. This is accomplished by 
using the Fox- Wolfram parameters [|26j. We require R2 = H2/H0 < 0.5. This cut has an efficiency of 99.38 ± 0.02% 
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while rejecting 29.0±0.2% of continuum events in the B decay kinematical region. The efficiency for this requirement 
is determined from our Monte Carlo simulation and agrees with the estimate derived from the data. 

Each charged track used to reconstruct a D or D* is required to be consistent with originating from the primary 
vertex. If the momentum of a track is greater than 0.250 GeV, we require that the z coordinate (along the beam) 
of the point of closest approach of the track to the beam-line be within 3 cm of the z coordinate of the event vertex 
and that the track's impact parameter with respect to the beam line be less than 5 mm. For tracks with momentum 
less than 0.250 GeV these requirements are loosened to 5 cm and 10 mm respectively. 

Particle identification requirements were imposed based on the specific ionization (dE/dx) measurements for the 
track, provided that more than 10 hits were recorded in the main drift chamber. The observed dE/dx had to be 
within three standard deviations of that expected for the particle species considered. 

Each photon candidate shower is required to lie within the good barrel region (cos 6 < 0.707) of the crystal 
calorimeter and to have a minimum energy of 30 MeV. Photon candidates are also required to be well separated from 
the extrapolated position of all charged tracks, and the lateral shape of the shower should be consistent with that 
expected from an electromagnetic shower. 

Candidate 7r° mesons are reconstructed from pairs of photon candidates. If the effective mass of two photons is less 
than 2.58<r away from the expected ty° mass, the combination is accepted as a n candidate and then is kinematically 
fitted to the nominal 7r° mass. 



B. Common procedure 

The inclusive B decay spectra and branching fractions are obtained by the subtraction of the candidate mass 
distribution below BB threshold (scaled by the ratio of luminosities and of the e + e~ annihilation cross section) from 
the candidate mass distribution on the T(4S) resonance. To illustrate the effect of continuum subtraction we show 
in Fig. ^ the (K~ir + ) effective mass distribution of D° candidates in the sample taken at the T(4S) resonance, in 
the whole momentum interval allowed in B decay, and the corresponding distribution from the sample taken below 
BB threshold, scaled by a factor of 2.08, the combined ratio of luminosities and of the e + e~ (non bb) annihilation 
cross sections at the two energies. 

In order to perform such subtraction it is convenient to use a scaling variable. We use the scaling variable x, 
defined as x — p/p ma x where p max = 4.950 GeV/c, the momentum for a D° produced in the reaction e + e~ — » D°D 
at a center of mass energy of 10.58 GeV (p m ax ~ 4.920 GeV '/c for the continuum sample, at a center of mass energy 
of 10.520 GeV). 

The momenta of the charmed mesons are measured in the T(45*) rest frame rather than in the B rest frame. Since 
the T(4£>) mass (10.58 GeV) is slightly above the threshold for B meson pair production (10.56 GeV), the B mesons 
are not at rest. The B momentum ranges from about 265 MeV to 355 MeV (± one standard deviation). This motion 
smears the value of x relative to what it would be if the B were at rest. A Monte Carlo study shows that the smearing 
in the variable x varies from ±0.013 to a maximum of ±0.020. Taking this effect into account, the maximum value 
of x for B decay to D is 0.506, and 0.496 for decay to D* . 

1. Spectra and the B branching fractions 

In order to measure the spectra, we divided our sample of charmed particle candidates into 10 x bins (20 for the 
D° where the statistics are high). For each x bin we generated the effective mass distributions of the D candidates 
(already selected as candidate D* decay products in the case of D* , see below) from "on resonance" and "continuum" 
data and performed bin-by-bin the scaled continuum subtraction, obtaining the mass distributions of the candidates 
from B decay. We then fitted this mass distribution, in each x bin, to the sum of the D signal and the various 
backgrounds. The shape of the signal, its parametrization and the different backgrounds will be described later. 

We performed an identical analysis (except for the non-existent continuum subtraction) on the Monte Carlo 
simulated events, to find the detection efficiencies as a function of x. The inclusive B decay branching ratios are 
calculated from the integral of the appropriate D (or D* ) efficiency-corrected spectrum. 

2. D* -tagging and "background-free" D samples 

We have studied the shape of the D signal and its momentum dependence, in data and in Monte Carlo simulation, 
using D°s from D* + — > D°ir + decays and D + s from D* + — > D + ir° decays. The D° sample selected this way is of 
known flavor because of the charge of the D* + . It was also used to obtain the shape and momentum dependence of 
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the "switched-mass" D° background discussed later in Sect. [V A. We obtain a very low background sample of D°s 
from D* + by selecting events in a ±2a region around the peak of the mass difference Sm = m(K~Tif 7r^~) — m(K~ 7rJ~ ) 
distribution (a ~ 0.8 MeV). We obtain a low background sample of D + s using a similar selection on the Sm = 
m{K-TT + Tv + 7v°) -m(K- 7r + 7r + ) distribution (a ~ 1.0 MeV). The Sm distributions for D candidates in the D mass 
signal region in the three channels considered in this analysis are shown in Figs. ^Q^(data). Most of the residual 
background in the Z?*-tagged D mass distribution is eliminated by subtraction of the Sm sidebands. 

Due to the fact that the signal in 8m is very close to threshold, it is difficult to choose sidebands that are wide 
enough to balance the number of background events under the Sm peak and far enough from the peak. Our choices 
for signal and sideband regions are shown in Table |. The signal region has a half width, and the sidebands a total 
width, of approximately twice the Gaussian a of the Sm signal. The position of the signal peak is obtained from the 
data. 

Before subtraction, the number of events in the Sm sidebands must be scaled to the estimated number of background 
events in the peak region. Two methods were used to estimate the scaling factor. The first method is to fit the Sm 
distribution with a smooth threshold function background and a double-Gaussian signal and calculate the scale factor 
by integration of the background function. The second method is to obtain the scale factor by fitting the D invariant 
mass distribution in the Sm peak and side band regions with a two-Gaussian D signal function and a polynomial 
background. The scale factor is then given by the ratio of the background levels under the D peaks. Fig. ^| was 
obtained by the second method. The results of the two methods agree and the scale factors are close to 1.1. The same 
procedure was also used for the analysis of the Monte Carlo sample. The candidate D mass distributions obtained 
from data by this D*-tagging procedure are shown in Figs. [| and ^. 



C. Estimate of the detection efficiency 

The detection efficiencies for the char m mesons studied were obtained analyzing the dedicated Monte Carlo sim- 



ulation samples described in Sect. [I A with the same procedures and selections used to analyze data. However, 
measurement errors are slightly underestimated in the simulation. In practice the only relevant parameter is the 
width of the reconstructed mass. Consequently, we analyzed the data using the parameters extracted from the data 
themselves, and analyzed the Monte Carlo samples using the parameters extracted from the Monte Carlo samples. 



D. Calculation of the B decay branching fractions 

The differential B decay branching fractions are calculated bin-by-bin from the following equations: 

dB(B - DX) x B(D) = rfiVp (6) 

^ n BB x e 

dB(B -> D*X) x B(D") x B(D) = ^ (7) 

^ n BB x e 

where dNo (dNr>*) is the yield of D (D*) in that bin, N B g is the number of of e + e~ — > BB events produced and e is 
the ^-dependent detection efficiency. We use CLEO results | p7|j28po[ ] for the D* and D absolute branching fractions 

B{D* + -> D°ty + ) = (68.1 ± 1.0 ± 1.3)% (8) 

B(D* + ~* D + ty°) = (30.8 ± 0.4 ± 0.8)% (9) 

B(D*° -> D°tt°) = (63.6 ± 2.3 ± 3.3)% (10) 

B{D° -> K'ty + ) = (3.91 ± 0.08 ± 0.17)% (11) 

B(D + -> K--K+TT+) = (9.19 ± 0.6 ± 0.8)% (12) 

The errors shown in our results for the differential spectra are relative, i.e. bin-to-bin errors. They do not include the 
error from smoothing the fitting function parameters and the errors of the detection efficiency. The overall statistical 
error on the B decay branching fraction is, however, derived from the integral of the spectrum when no smoothing 
was performed. 
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E. Systematic Errors 



In this section we discuss the sytematic uncertainties of the B — > D X and B — > D* X analyses. Uncertainties 
related to specific decay channels are deferred to the relevant sections. We estimate a 1% uncertainty in charged 
track detection efficiency (2% for p(n) below about 200 MeV/c) and 5% in 7r° detection p3|. 

These errors are coherent for each kind of track and hence should be multiplied by the respective number of tracks 
in the decay under study. However, we obtain the B decay branching fractions by dividing the product of branching 
fractions (Eqs. ^ and (Q)) by the D (and D*) decay branching fractions measured in our own experiment. Hence, to 
the extent that the kinematical configurations and the data set used overlap, most tracking errors are eliminated. We 
estimate a residual systematic error of 0.5% per charged track (1.0% if below 200 MeV/c) and 2% for tt° whenever 
there is compensation. 

We have studied the effect on the branching fraction of the track quality, geometry cuts and event shape cuts by 
succesively removing them and then measuring the effect on the B branching fractions. The errors specific to each 
final state will be quoted in the appropriate sections. 

In measuring B(B — > D* X) the subtraction of the background due to association of a true D and a random it is 
performed by subtracting the D candidate mass distribution of the Sm sidebands from that in the Sm peak region. 
The result is slightly dependent on the choice of the width of the Sm intervals (Table |). We varied the widths of 
these intervals by ±0.5 MeV in each case and derived a relative uncertainty in the B branching fraction of 0.9% for 
B -> D*+ and of 1.2% for B -> D*°. 

In the case of B — > D* decays, the branching fraction is also sensitive to the value of the scale factor used in the Sm 
sideband subtraction This effect contributes a systematic error which is estimated to be the change in the branching 
fraction corresponding to a la change in the scale factor. 

The value of N B g in Eqs. ^ and (^) is affected by a relative error of 1.8% (Sec. |n|). Finally, we take into account 
the statistical error on the efficiency, calculated from the Monte Carlo simulation. All the above errors are combined 
in quadrature to give the relative systematic error on our measurements. 

The third error in the B branching fractions arises from the propagation of errors in the D and D* decay branching 
fractions, in so far as they are not coherent with the errors in this analysis. 



IV. INCLUSIVE B -> D° DECAY 

The continuum-subtracted m{K~ distribution in any given D° momentum interval is the sum of the D° signal 
plus various backgrounds. This is illustrated in Figs. ^ and |^ which show the m(K~n + ) distributions for one of 
the lowest momentum intervals (0.050 < x < 0.075), and for a higher momentum interval (0.325 < x < 0.350), 
respectively. We see a prominent signal and three backgrounds: (i) a combinatorial background; (ii) the "switched 
mass" background that contributes to the signal region and (iii) two "satellite bumps" on either side of the signal. 
Backgrounds (ii) and (iii) complicate the fit to the signal region and its immediate vicinity. It is not possible to 
establish the shape of the combinatorial background without a reliable knowledge of the shapes and amounts of 
backgrounds (ii) and (iii) . A detailed discussion of the signal shape and of these backgrounds is presented below. 



A. Switched mass background 

D° —> n~K + events may be misidentified as D —* K~n + when there is insufficient K/tt discrimination. This 
background complicates the extraction of the D° yield, particularly at low D° momenta (p < 700 MeV) where this 
"switched mass" background distribution is similiar in shape to the signal and peaks under it (see Fig. Q and ^|). 
For a D° at rest and without particle identification the area of this background within 2a of the signal peak would 
equal the signal area. However, this background, within such limits, quickly decreases to about 10% of the area of 
the signal for D° momentum above 1 GeV and decreases further at higher momenta. For xd > 0.50, the invariant 
mass distribution of these doubly misidentified events becomes so broad that it can be absorbed in the polynomial 
parameterization of the combinatorial background. The amount of this type of background is strongly affected by 
the use of particle identification. 

We studied the shape and relative size of this background and its dependence on D° momentum in the data and in 



the Monte Carlo simulation using D* + -tagged D°s (Sec. IIIB 2), whose flavor is tagged by the charge of the tv + . We 



produced "switched mass" distributions, by interchanging the mass assignments of the K and it tracks for each event, 
provided that the inverted assignements were consistent with our particle identification criteria. We then subtracted 
the "switched mass" distribution in the Sm sidebands from that in the Sm peak region. Figs. ^ and ^ show the 
resulting m(ir~ K + ) distributions for a low D° momentum interval, and for a higher D° momentum interval. We also 
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calculated, as a function of momentum, the area of this background relative to the area of the correctly reconstructed 



IIIB2 



D signal also obtained as described in Sec. 

With our statistics it is difficult to smoothly parameterize the shape of this background as a function of momen- 
tum. Figs. ^ and ^ show, however, that the switched mass distributions obtained with the MC-tag analysis of the 
gener ic Mon te Carlo simulation provide an excellent, parameter-free description of the D* + -tagged data distributions 
(Sec. 
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The preceding analysis is performed both on the data and on the dedicated Monte Carlo samples. The momentum 
dependence of the ratio of the switched mass distribution area relative to the area of the signal is shown in Fig. jTl] 
as obtained from the MC-tag sample. Similar dependences are obtained using the D* + -tagged samples of data and 
dedicated Monte Carlo. We have checked that changing from one dependence to another has minimal impact on the 
final result. 



B. Background from misidentified D° -> K~K+, D° -> tt~tv + and D° -> K'tt + tt' 



The Cabibbo-suppressed D° decay modes, D° — > K~ K + and D° — > iv~-k + , produce asymmetric enhancements on 
opposite sides of the signal region when one of the tracks is misidentified as a ty or K, respectively (see Fig. |^ and ^|). 
The K~iv + effective mass distribution from the decay D° — ► K~tv + tt° introduces an additional small background at 
the lower edge of our mass spectrum. Switching the kaon and pion masses in these decays produces an additional 
background at masses well below the D° mass but with tail up to the D° mass. If these contributions are not 
included in the shape of the background, it is possible to overestimate the level of combinatorial background and 
thereby underestimate the signal. Not taking into account the presence of these enhancements both in fitting the 
data and in fitting the dedicated Monte Carlo for obtaining the efficiency, results in a change of the order of 1% in 
the branching fraction. Also in this case we have found that the momentum dependence of the shape and of the area 
ratio obtained by the MC-tag procedure from the generic Monte Carlo is quite good and superior to that obtained 
using complicated analytical parametrizations. The size of these backgrounds relative to the signal was taken from 
the Monte Carlo simulation, but we have checked that letting their normalization float in the fit did not appreciably 
alter the result. 



C. Combinatorial background 



The combinatorial background, which is nearly flat in the m(K ir + ) region of interest, is parameterized by a 
second order polynomial. 



D. Raw D spectrum and the shape of the signal 



It is very important to obtain an accurate representation of the signal shape and specifically of its tails, because 
of the presence of the mass-switched background and its correlation with the signal. We used two diferent signal 
parameterizations. 



1. Analytical parameterization 



A Gaussian function does not give a sufficiently accurate parameterization of the D° signal. Track measuring 
errors are different in different events because of the geometrical orientation of the D° decay products in the detector 
and their overlap with other tracks. In order to take into account this variation of the errors, one possibility is to use 
a double Gaussian with the two Gaussians constrained to have the same mean. The parameters of the signal shape 
are the mean, the width of the narrower Gaussian, ax, the ratio of the widths of the wider to the narrower Gaussian, 
ail ax, and the ratio of the area of the wider Gaussian to the total area, Az/Atot- There are strong correlations 
among these three shape parameters. It is then difficult to obtain a smooth D° momentum dependence for them if 
they are allowed to float when fitting the signal in each momentum bin. 

We obtained satisfactory results for these signal shape parameters by (i) combining "on resonance" and "contin- 
uum" D° candidates together to obtain hig her sta tistics, especially close to the kinematic limit, or (ii) by using the 
D*-tagged, background-free D° signal (Sec. Ill B 2 ). We chose the set that gave the best fit to the mass distributions 
and smaller errors for the parameters. In order to minimize the bin-to-bin statistical fluctuations, we re-binned the 
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D° sample in only 9 momentum intervals (0.025 < x < 0.475) and studied the position and shape of the D° signal 
peak as a function of momentum. We smoothed the momentum dependence of the signal shape parameters so ob- 
tained with polynomial functions of the momentum. We then obtained the raw D° spectrum by fitting the m(K~w + ) 
distribution in 20 x bins between 0.0 and 0.5 keeping a^/ox and Ai/Atot fixed at the smoothed values. o\ is allowed 
to float and varies from about 6 MeV at low momentum to about 8 MeV at the maximum momentum; oijox and 
A^IAtot are approximately constant, the first at about 2.7 and the second at about 0.15. 

The same procedure was followed to analyze the events from the Monte Carlo. The comparison of the data and 
Monte Carlo shape parameters shows that the Monte Carlo simulation underestimates the track measuring errors 
(reflected essentially in o\) by about 10%. 



2. Use of Monte Carlo simulation with track tagging 

An alternative way of parameterizing the D° signal shape and a check of the double-Gaussian parameterization, is 
to us e the histograms provided by our generic Monte Carlo simulation and tagging the reconstructed tracks (MC-tag, 



Sect. II A). In constructing these histograms we have used the same selection criteria used in the analysis of data 
and of dedicated Monte Carlo already described. 

We have just seen that our Monte Carlo simulation underestimates by about 10% the overall width a of the D° 
signal over the whole D° spectrum. It is however likely that the Monte Carlo simulation reproduces more accurately 
the dependence of a on track position, orientation and overlap with other tracks, all factors that cause the non- 
Gaussian spread in measuring errors seen in the data. We then expect that the shape of the D° signal obtained 
by tagging the Monte Carlo tracks will be a good representation of the data if corrected for the overall width. We 
have expressed this correction by one parameter (then fitted on the data) that changes the width of the signal mass 
distributions without otherwise altering its shape. 

We used the histograms obtai ned by this procedure as fitting functions of the "background-free" signal distributions 



obtained by D* -tagging (Sect. IIIB 2). We obtained excellent fits, altough not superior to the double- Gaussian fits. 
We obtained D° spectra using these distributions as fitting functions. The D° yields differed from those obtained 
with the double-Gaussian by less than a tenth of the statistical error. 



E. Results 



Particle identification (Sect. [II A) reduces the combinatorial background by about a factor of 2 and improves the 
statistical error while slightly increasing the systematic error due to the small uncertainty on the accuracy of the 
simulation of the ionization measurements. We performed our analysis both with and without particle identification 
requirements and obtained results that were in good agreement. The results presented here were obtained using 
particle identification and the double-Gaussian parameterization of the signal. In Figs. [t| and|l| we show the fitted 
mass distributions in the 20 x bins. Summing the spectrum over the 20 bins for 0.0 < xd < 0.5 gave a raw yield of 
62, 648 ± 1, 394 D° from B decay. 

The dedicated Monte Carlo simulated events were analyzed as real events to produce the histogram of detection 
efficiency versus momentum that was smoothed and used to correct the raw spectra (Fig. ^). The detection efficiency 
is nearly independent of momentum in the region of interest, except at the lowest momenta where there is an increase 
of angular acceptance due to the near alignment of the two D° decay products. 

The continuum-subtracted, efficiency-corrected, inclusive D° momentum spectrum in B — > D°X decay is shown 
in Fig. |l5| Summing this spectrum over the interval < x < 0.5, we obtain the efficiency corrected yield of 
108, 507 ± 2, 407 D°s (resulting in a spectrum averaged detection efficiency = 0.578) out of 4.3314 million B decays. 
This gives the product of branching fractions, 

B b Bd = (2.51 ± 0.06)%, (13) 

where B B = B(B -> D°X) and B D = B(D° -> K~-k + ) + B(D° -> K + n~) and the error is only statistical. 

Using the CLEO results BjD° -> R-TT+) = (3.91 ± 0.08 ± 0.17)% g] and F{D° -> K + -r-)/Y{D° -> K~tt + ) = 
(0.0077 ± 0.0025 ± 0.0025) wk, we obtain from Eq. ([u| the branching fraction 



x Our procedure finds both the Cabibbo-allowed, D° — > K tv + , and the doubly Cabibbo-suppressed, D° — » K + ir , 
decay modes, so we must divide by the sum of the two decay branching fractions. 
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B(B -> D°X) = 0.636 ± 0.014 ± 0.019 ± 0.018. 



(14) 



The first error is statistical, while the second is the systematic error. The contributions to the second (systematic) 
error, as percentages of the branching fraction are given in table [□]. The third error is due to the error on the absolute 
branching fraction B(D° — > K~tt + ). 

As part of a different project |B3[, an analysis of the B — > D°X decay was carried out over a somewhat increased 
data sample and with more stringent particle identification requirements, using also time of flight information. If 
both the D° — > K~ir + and ir~ K + hypotheses were acceptable, but one gave a particle ID Xworse > 4 + xLtten the 
hypothesis with the worse x 2 was rejected. This selection in fact reduced the switched mass background to zero. The 
result of this analysis for B(B — > D°X) is fully consistent with our result. 



V. INCLUSIVE B -> D+ DECAY 



The complications from satellite peaks or switched mass backgrounds are not present in the measurement of the 
D + spectrum and other kinematic reflections are quite small (they will be discussed in the next subsection). It is then 
possible to estimate the shape of the combinatorial background by fitting the K~tv + -k + mass distribution excluding 
the signal region. However, it is difficult to determine the momentum dependence of the signal shape parameters 
given the large combinatorial background and low signal to background ratio. The K~tt + tv + mass distributions 
summed over all momenta of interest for the data at the T(45*) and those below BB threshold (scaled according 
to luminosity and square of center-of-mass energy ratios) are shown in Fig. |l(| If we parameterize the signal as a 
double Gaussian, a modest change in the curvature of the combinatorial background in the m(K~ir + ir + ) distribution 
strongly correlates with a wider and larger second component of the double Gaussian signal. 

We analysed the m(K~ir + ir + ) distributions using the double Gaussian parameterization of the signal. As in the 
analysis of the D° spectra, we used D + tagged by the D* + — > D + ir° decay to determine the momentum dependence 
of the signal parameters. We also performed the analysis using the simple Gaussian signal parameterization. This is 
quite adequate to fit the data, given our statistical errors. We found that the results from the two different signal 
parameterizations are the same within a small fraction of the statistical error. Here we shall present the results 
obtained with the simple Gaussian parameterization. 



A. Background from Df Decays 



We have small backgrounds from Df decays where a K + is misidentified as a 7v + . We took this into account using 
the following procedure. We generated BB Monte Carlo events that contain at least one Df, which decays into 
K~ir + K + directly or through the K* K + or 4>tt + resonance channel. We processed them through the simulation 
of our detector and track reconstruction. We selected events that passed the selection criteria for the decay D + — > 
K~it + it + . We plotted these (fake) m( K~-k + -k + ) distributions for each momentum bin and normalized them using 
our knowledge of the differential B(B — » Df X) [B4J and of the D s decay branching ratios. We added these histograms 
to the combinatorial background and signal functions when fitting the data. Taking into account this feed-down we 
changed our result for B(B —+ D + X) by AB = —0.0021 ± 0.0002. We have examined other possible feed-downs and 
verified that they do not affect our measurement. 



B. D spectrum 



We fitted the continuum subtracted sample of D + — > K~tt + iy + candidates. We divided the sample in 20 momentum 
bins between 0.00 < x < 1.00 and performed the bin-by-bin fit of the m(K~ -k + -k + ) distributions using the signal and 
background functions described above. In Fig. ^| we show the fitted mass distributions in the 10 x bins. The sum 
over the 10 bins for 0.0 < xo < 0. 5 g ave a raw yield of 35, 804 ± 1, 297 D + from B decay. The detection efficiency 
as a function of x is shown in Fig. |l£|. We have smoothed the x dependence in different ways. The resulting average 
efficiencies differ from one another by 0.18% at most. The error introduced by the choice of the smoothing function 
has been taken into account in the systematic error. The smoothing actually used is shown in Fig. O. 
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C. Result 



The efficiency corrected D + spectrum is shown in Fig. Summing this spectrum over the interval 0.0 < xr> < 0.5, 
we obtained an efficiency corrected yield of 93, 493 ± 3, 386 D + (resulting in a spectrum averaged detection efficiency 
= 0.383) that gives the product of branching fractions, 

B b Bd = (2.16 ± 0.08)%, (15) 

where Bb = B(B — > D + X), Bd = B(D + — > K~tt + -k + ) and the error is statistical only. 

Using CLEO's result for the absolute branching fraction B{D+ -> K~n + n + ) = (9.19 ± 0.6 ± 0.8)% @ gives: 

B{B -> D + X) = 0.235 ± 0.009 ± 0.009 ± 0.024. (16) 

The first error is statistical, the second is systematic, and the third error is due to the statistical error on B(D + — > 



K - 7r 7T ) and to that part of its systematic error that propagates to our measurement. Table [II gives the components 
of the second error as percentages of the branching fraction. 

Also this analysis, as in the B — > D° X case, has been carried out with different procedures and selection criteria 
that produced results within a fraction of the systematic error of the one reported here. 

We will show later that the B — > D* + X branching fractions measured by reconstructing D* + through the decay 
chains D* + — » D°7v + — > (K~ tv + )tt + and D* + — » D + 7v° — » (_ftT~7r + 7r + )7r° are consistent with each other, their ratio 
being 1.01 ± 0.09. This is a consistency check of the results given above. 



VI. INCLUSIVE B -> D* DECAY 



We selected D* candidates by combining a 7r of the appropriate charge with a D candidate and then imposing 
requirements on 8m, the difference between the mass of the D* candidate and that of the D candidate. However, this 
sample also contains background composed of pions randomly associated with a correctly reconstructed D meson . 
We eliminated this "fake" D* background by subtracting properly scaled 5m "sidebands" as described in Sec. [II B 2] , 
The mass distribution of the D candidates so obtained was fitted to find the D yieid as in the previous sections. 



A. Inc[usive B — > D Decay 



We se[ected D° candidates in the decay mode D° — > K~ir + (as described in Sect. fV\ ) and combined them with a 
7T° to form D*° candidates. The D*° yieid in each momentum bin was found by fitting the m(K~ tt + ) distribution 
obtained with the 5m seiection and sideband subtraction procedure. This decay channei is affected by the same 
backgrounds as the D° channei. To fit the m(K~ mass distributions, we used the parameterization and procedures 
from the analysis of B — > D° + X described in Sect. However, the large number of free parameters involved with 
the multi-Gaussian parameterization was not suitable with the lower statistics available in this channel. We have 
thus used the histograms of the satellite peaks and of the switched K — n distributions from Monte Carlo tagging 
while for the signal we used either a single Gaussian or, the histogram from Monte Carlo tagging. Both procedures 
gave B(B — > D*° + X) values differing by a small fraction of the statistical error. We report here the results from 
the single Gaussian parameterization of the signal. In Fig. ^ we show the fitted mass distributions in the 10 x 
bins. The sum over the 10 bins gave a raw yield of 3, 539 ± 175 D*° from B decay. After correcting bin-by-bin for 
detection efficiency (Fig. ^), we obtained the inclusive D*° spectrum shown in Fig. ^3| Fig. ^ shows the continuum 
subtracted K~tt + effective mass distribution of D*° candidates from the decay chain B — > D^X — > (D°ir°)X in the 
momentum interval (0.0 < xd' < 0.50). 

As in the previous cases, by summing the spectrum in the interval 0.0 < x < 0.5 we found a corrected yield of 
26, 840 ± 1, 331 D*° from B decay (resulting in a spectrum averaged detection efficiency = 0.132) and the product of 
branching fractions, 

B b B* d Bd = (0.620 ± 0.031)% (17) 

where B B = B(B -> D*°X), B* D = B(D*° -> D°ty°), Bd = B(D° -> K--K+) and only the statistical error is shown. 

Using the CLEO results B{D*° -> D°ty°) = (63.6 ± 2.3 ± 3.3)% @ and B{D° -> K~-k+) = (3.91 ± 0.08 ± 0.17)% 
||, and F(D° -> K + -k~)/T(D° -> K~n + ) = (0.0077 ± 0.0025 ± 0.0025) Q, from Eq. @ we obtain the branching 
fraction: 

B(B -> D*°X) = 0.247 ± 0.012 ± 0.018 ± 0.018. (18) 
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The first error is statistical. The second error is systematic and includes the components listed in Table WWl These 
are given as percentages of the branching fraction itself. The third error is determined from the statistical error on 
B{D° -» K--K+). 



B. Inclusive B — > D* + Decay 

We have analysed the inclusive decay B — > D* + in two D* + decay modes. The D* + — > D°tt + channel has 
significantly more events because of the larger D* branching fraction and higher detection efficiency. However, 
the detection efficiency is a steep function of the D* + momentum (Fig. because of the short range of the low 
momentum 7v + and absorption in the beampipe. The detection efficiency for the charged pion is nearly zero for 
xd* < 0.15. Summing over the spectrum gives B(B — » D* + X) in the chargedpion mode for xd* > 0.15. The 
D* + — > D + 7T° channel, which has fewer events and much larger backgrounds (Fig. fa) has an efficiency which is nearly 
constant with momentum (Fig. ^) and is the only source of information for the low momentum region. We will 
separately describe the analyses using the two D* + decay modes and then discuss how to combine them to obtain 
B{B -> D*+X). 



1. Using the D* + — » D°ir + decay channel 

Candidate D° mesons were reconstructed in the D° — > K~ty + decay mode (as described in Sec. f\/\ ) and combined 
with a 7r + (referred to as the "slow ir") to form D* + candidates. In order to maximize the detection efficiency of 
the slow 7r + we did not require dE/dx information to be available for this track. The D* + yield in each of the 
7 momentum bins (0.15 < xd* < 0.50) was found by fitting th e m(K ~ ir + ) distribution obtained through the 5m 



selection and sideband subtraction procedure described in Sec. [ill B 4 In fitting the D° peak, we used the same 
parameterization that was used in the analysis of B — > D° + X, but no background due to double misidentification 
of D° -> R- 7r + is present in this mode. In Fig. Ea we show the fitted mass distributions in the 7 x bins. Summing 
the spectrum in the interval 0.15 < xd* < 0.5, gave a raw yield of 4, 214 ± 120 D* + . 

Correcting bin- by-bin for detection efficiency (Fig. |24|), we obtained the inclusive D* + spectrum from B decay 
shown in Fig. ^(j. The efficiency corrected yield for xd* > 0.15 is 24, 691 ± 902 D* + from B decay (resulting in a 
spectrum averaged detection efficiency = 0.171). This gives the product of branching fractions, 

BbB* d Bd{xd* > 0.15) = (0.570 ± 0.021)% (19) 

where B B = B(B -> D*+X), B D = B{D*+ -> £>%+), B D = (B(D° -> R--K+) + B(D° -» K + iz~)) and the error 
quoted is only statistical. 

Fig. bjj shows the continuum subtracted K~tv + effective mass distribution of D* + candidates from the decay chain 
B — > X — * (D°tt + )X in the momentum interval interval (0.15 < xd* < 0.50), after 5m sideband subtraction. 



2. Using the D* + — > D + -k° decay channel 

Candidate D + s are reconstructed in the decay mode D + — > K~tt + tt + (as described in Sec. |y]) and combined with 
a 7T° to form D* + candidates. The D* + yield in each of the 20 momentum bins (0.0 < xd* < 1.0) was found by 
fitting t he m(K ^-K + 7v + ) distribution obtained from the 5m selection and sideband subtraction procedure described 



in Sec. [II B 2 



The D peak was fitted using the same parameterization and procedure used in the analysis of 
B — » D + X. Summing the spectrum over the interval 0.0 < xd* < 0.5 gave a raw yeld of 2,925 ± 250 D* + . After 
correcting for detection efficiency (Fig. |28[ ) and summing over x bins, we obtained the inclusive D* + spectrum from 
B decay shown in Fig. |2i| The corrected yield is 27, 683 ± 2, 339 D* + from B decay (resulting in a spectrum averaged 
detection efficiency = 0.106) and the product of branching fractions is, 

BbB d Bd = (0.639 ± 0.054)% (20) 

where B B = B{B -> D* + X), B* D = B{D* + -> D + n°), B D = B{D + -> K~n + n + ) and the error quoted is statistical 
only. 

Fig. ^ shows the continuum subtracted K~7v + tv + invariant mass distribution of D* + candidates from the decay 
chain B — > D* + + X — > (D + w°)X in the whole momentum interval interval (0.0 < xd* < 0.50). 
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C. Combined results for the inclusive B — > D* + Decay 

Using the corrected differential branching fractions dB{B — > D* + X)/dx obtained in two independent D* decay 
modes, we combined them as follows. In the momentum region 0.0 < xd* < 0.15 we use the only measurement 
available, that from the D* + -> D + ty° decay mode, B(B -» £>*+X)(0.0 < x D * < 0.15) = 0.031 ± 0.009 ± 0.0025 ± 
0.0027. In the momentum region 0.15 < xd* < 0.50 we calculated the weighted average of the two measurements. 
The resulting spectrum is shown in Fig. 

The sum over all x bins gives the branching fraction, 

B(B -> D* + X) = 0.239 ± 0.011 ± 0.014 ± 0.009. (21) 

The first error is statistical and is dominated by the error on the D* branching fraction for 0.0 < id* < 0.15. The 
second error is systematic and includes the components listed in Table [y]. These are quoted as percentages of the 
branching fraction itself. The third error is propagated from the statistical error on the D°, D + and D* + decay 
branching fractions. 



VII. POLARIZATION 



The polarization of D* mesons has been predicted for semileptonic B decays and for two body hadronic decays 
36) . A model dependent estimate of momentum dependence of the polarization for directly produced D* mesons is 
available for inclusive decays ||. CLEO has previously measured the D* + polarization [^B). Here we present a new 
measurement of the D* + polarization with higher statistics. However, the present detector is operated with a higher 
magnetic field that makes impossible to extend the measurement to low D* + momenta for the D* + — » D°ir + decay 
mode. We also present the first measurement of D*° polarization in inclusive B decays. In this case, the measurement 
can be extended to the lowest momenta without difficulty. These polarization measurements served also as a check 
of the accuracy of our Monte Carlo simulation of B decay. 

The polarization as a function of x is determined from the distribution of the D* decay angle, 9. This is the angle 
between the direction of flight of the D* in the laboratory frame and the direction of the daughter D° in the D* rest 
frame. The distribution of this decay angle can be parameterized as 

dW 3 , 

-t— — — - = 1 + a cos 2 9} (22) 

dcos9 6 + 2a L J y ' 

This is equivalent to the expression in terms of the spin-density matrix element poo: 

W (cos 6») = | [(1 - poo) + (3poo - 1) cos 2 9] (23) 

Here, the spin-density matrix is determined in a coordinate system with the quantization axis along the direction of 
motion of D* in the laboratory frame. The element poo is the probability for the D* to be in the J z = state. 
The polarization parameter a is related to the longitudinal and transverse decay rates of D* as 

Tl i 

When a is close to its lower bound of — 1, the distribution is sin 2 9 implying transverse polarization while large 
values of a imply longitudinal polarization. 

Only the D* — > D°n, D° — > K~tv + decay mode is used. We applied the same selection criteria that were used in the 
branching fraction measurements to obtain the D* samples. The sample was divided into 5 intervals in 0.0 < x < 0.5 
and in 5 intervals in cos 9. Each of the Kn distributions was fitted to a Gaussian shape plus polynomial background. 
We repeat the same analysis procedure on the Monte Carlo simulated data to find the dependence of efficiencies on 
x and cos 9. 

Figures ^ and |3^ show the efficiency corrected and background subtracted cos 9 distributions for D* + and D*° 
respectively. We have compared these distributions with those predicted by our Monte Carlo simulation. The 
simulation appears to correctly model the data. This checks our results on spectra and B decay branching ratios 
because an incorrect simulation of the polarization may result in an incorrect determination of the detection efficiency. 

Figure ^ shows a as a function of the scaled momentum variable x for D* + and D*° . 
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VIII. DISCUSSION AND SUMMARY 



In Table [Vl| we compare our results with those from previous measurements. In order to make the comparison 
independent of D* and D decay branching fractions used in the different experiments, we give the product of the 
branching fractions. 

Only the statistical error is reported here for our current results (third column) because the systematic error on 
these products of branching ratios cannot take advantage of some cancellations when we divide them by D and D* 
branching ratios measured in our own experimant. Our results are consistent with previous measurements, except 
for the B -» D°X case where our branching fraction is significantly higher. The branching fraction B(B — > D*°X) is 
measured here for the first time. The mode B(B D* + X) is reconstructed in the D* + — » D + tt° decay mode for the 
first time. The result for B(B — > D* + X) in the D* + — > D°tt + decay mode refers only to xd > 0.15 and the previous 



one from CLEO only to xd > 0.10. Table VII summarizes the results on the inclusive B decay branching fractions. 
From the measurement of B(B — > D*°X) we can determine the ratio 

B(B -> D*°X) 

b\b^d* + x) =im± a07 ± a09 ± a08 

which is consistent with the naive expectation of 1.00. 



Using these new measurements, and previous ones shown in Table VIII , we can now calculate the average number 
of charm quarks produced in B decay, 

(n c ) = 1.10 ±0.05, 

This value is consistent with the naive expectation of 1.15, but it is considerably lower than the value (~ 1.30) 
required to account for the low value of the B semileptonic branching ratio in models where the channel b — > ccs is 
enhanced §00000- 

ACKNOWLEDGEMENTS 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity and running 
conditions. J.P.A., J.R.P., and I.P.J.S. thank the NYI program of the NSF, M.S. thanks the PFF program of the 
NSF, G.E. thanks the Heisenberg Foundation, K.K.G., M.S., H.N.N. , T.S., and H.Y. thank the OJI program of 
DOE, J.R.P., K.H., M.S. and V.S. thank the A.P. Sloan Foundation, and A.W. and R.W. thank the Alexander von 
Humboldt Stiftung for support. M.S. is supported as a Scholar of the Cottrell Research Corporation. This work 
was supported by the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences and 
Engineering Research Council of Canada. 



REFERENCES 

[1] For a review of the experimental and theoretical status of B meson decays, see"B Decays", S. Stone ed., 2nd Ed., 
World Scientific, Singapore, 1994; "B Mesons", T.E. Browder and K. Honscheid, to appear in Progress in Nuclear 
and Particle Physics, Vol. 35, 81-220 (1995). 

[2] A.F. Falk, M.B. Wise, I.Dunietz, Phys. Rev. D 51, 1183 (1985). 

[3] M. Wirbel and Y.-L. Wu Phys. Lett. B 228, 430 (1989). 

[4] I.I. Bigi, M. Shifman, N.G. Uraltsev and A.I. Vainshtein, Phys. Rev. Lett. 71, 496 (1993); B. Blok, L. Koyrakh, 

M. Shifman and A.I. Vainshtein, Phys. Rev. D49, 3356 (1994). 
[5] A. Manohar and M.B. Wise, Phys. Rev. D49, 1310 (1994). 
[6] T. Mannel, Nucl. Phys. B413, 396 (1994). 

[7] A.F. Falk, Z. Ligeti, M. Neubert and Y. Nir, Phys. Lett. B326, 145 (1994); L. Koyrakh, Phys. Rev. D49, 3379 

(1994); S. Balk, J.G. Korner, D. Pirjol and K. Schilcher, Z.Phys. C64, 37 (1994). 
[8] I.I. Bigi, N.G. Uraltsev and A.I. Vainshtein, Phys. Lett. B293, 430 (1992); Erratum, Phys. Lett. B297, 477 (1993); 

I.I. Bigi, B. Blok, M. Shifman, N.G. Uraltsev and A.I. Vainshtein, 7th Meeting of the Division of Particles Fields 

of the APS (DPF 92), Batavia, IL, 10-14 Nov 1992, published in DPF Conf.l992:610-613. 
[9] B. Blok and M. Shifman, Nucl. Phys. B399, 441 (1993); B. Blok and M. Shifman, Nucl. Phys. B399, 459 (1993); 

N. Bilic, B. Guberina and J. Trampetic, Nucl. Phys. B248, 261 (1984); M. Voloshin and M. Shifman, Sov. J. Nucl. 

Phys. 41, 120 (1985); M. Voloshin and M. Shifman, Sov. Phys.-JETP 64, 698 (1986); V.A. Khoze, M. Shifman, 

N.G. Uraltsev and M. Voloshin, Sov. J. Nucl. Phys. 46, 112 (1987). 



15 



E. Bagan et al, Nucl. Phys. B 432, 3 (1995); Phys. Lett. B 342, 362 (1995); Phys. Lett. B 351, 546 (1995), [E: 
Phys. Lett. B 374, 363 (1996)]. 

M. Neubert and C.T. Sachrajda, "Spectator Effects in Inclusive Decays of Beauty Hadrons", preprint CERN- 



TH/96-19 ( fiep-ph/9603202D , submitted to Nucl. Phys. B. 

G. Buchalla, I. Dunietz and H. Yamamoto, Phys. Lett. B364, 188 (1995). 

B. Blok and T. Mannel, Mod. Phys. Lett. All, 1263, (1996). 

I.L. Grach, I.M. Narodetskii, S. Simula and K.A. Ter-Martirosyan, "Exclusive and Inclusive Weak Decays of the 



B-Meson", Istituto Superiore di Sanita', preprint no. INFN-ISS 95/17 (tiep-ph/9603239) submitted to Nucl. Phys. 
B. 

D. Bortoletto et al. (CLEO), Phys. Rev. D 45, 21 (1992). 

H. Albrecht et al. (ARGUS), Z. Phys. C 52, 353 (1991). 

B. Barish et al. (CLEO) Phys. Rev. D 51, 1014 (1995). 

Y. Kubota et al. (CLEO), Nucl. Inst, and Meth. A 320, 66 (1992). 

D. G. Cassel et al, Nucl. Inst, and Meth. A 252, 325 (1986). 

E. Blucher et al, Nucl. Inst, and Meth. A 249, 201 (1989). 
R. Morrison et al. (CLEO), Phys. Rev. Lett. 67, 1696 (1991). 
D. Bortoletto et al, Nucl. Inst, and Meth. A 320, 114 (1992). 

T. Sjostrand, Comp. Phys. Comm. 39, 347 (1986), T. Sjostrand, M. Bengston, Comp. Phys. Comm. 43, 367 
(1987). 

Our Monte Carlo simulation does not include the possibility of B — > DDKX decay. Since our detection efficiency 

has been obtained bin-by-bin in the D momentum, this omission does not affect our result. 

R. Brun et al., GEANT 3.14, CERN DD/EE/84-1. 

G.C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

F. Butler et al. (CLEO), Phys. Rev. Lett. 69, 2041 (1992). 

D. Akerib et al. (CLEO), Phys. Rev. Lett. 71, 3070 (1993). We use the value 3.91% before radiative corrections 
because our D° yield is not corrected for radiative corrections. 
D. Cinabro et al. (CLEO), Phys. Rev. Lett. 72, 1406 (1994) 
R. Balest et al. (CLEO), Phys. Rev. Lett. 72, 2328 (1994). 
Review of Particle Properties, Phys. Rev. D 54, 1 (1996). 

For a discussion of charged track and 7r° detection efficiency and their errors see ref. |l^] 
CLEO work in preparation. 

D. Gibaut et al. (CLEO), Phys. Rev. D 53, 4734 (1996). 

CLEO internal note CBX 96-43 

J.L. Rosner Phys. Rev. D 42, 3732 (1990). 

M. Zoeller, Ph.D. dissertation, S.U.N.Y. at Albany, (1994) 

D. Cinabro et al. (CLEO), paper contributed to the ICHEP conference, Glasgow, 1994, preprint CLEO CONF 
94-08, and J.P. Alexander et al. (CLEO) Phys. Rev. Lett. 74, 3113 (1995), E: Phys. Rev. Lett. 75, 4155 (1995). 
M.B. Voloshin, Phys. Lett. B385, 369 (1996). 
R. Balest et al. (CLEO) Phys. Rev. D52, 2661 (1995). 

TABLES 



TABLE I. 8m cut and width of the sidebands (in MeV). 



channel 


width 


signal region 


lower-band 


upper-band 




±1.5 


144.1-147.1 


141.1-142.6 


148.6-150.1 


D'+ -> D+7Y° 


±1.5 


139.1-142.1 


136.1-137.6 


143.6-145.1 




±2.0 


140.6-144.6 


137.6-139.6 


145.6-147.6 



1G 



TABLE II. Relative Systematic Errors on B(B -> D°X) 



residual track finding efficiency uncertainty 1.0% 

variation of track quality and geometrical cuts 0.5% 

choice of signal shape parameterization 2.0% 

error in particle identification efficiency 0.8% 

error in number of B and B 1.8% 

statistical error on efficiency 0.4% 

Total T0% 



TABLE III. Relative Systematic Errors on B(B -> D+X) 



residual track finding efficiency uncertainty 1.5% 

variation of track quality and geometrical cuts 1.2% 

choice of signal shape parameterization 2.0% 

error in number of B and B 1.8% 

statistical error on efficiency 0.7% 

uncertainty in background shape 1.7% 

smoothing of the efficiency vs x 0.14% 

estimate of kinematical reflections 0.10% 

Total 3.8% 



TABLE IV. Relative Systematic Errors on B{B -> D*°X) 



residual particle finding efficiency uncertainty 5.1% 

variation of track quality and geometrical cuts 0.5% 

choice of signal shape parameterization 2.0% 

error in number of B and B 1.8% 

statistical error on efficiency 2.1% 

choice of the 5m region widths 1.2% 

scale factor in the 8m sideband subtraction 3.3% 

Total TWo 



TABLE V. Relative Systematic Errors on B(B -> D* + X) 



residual track finding efficiency uncertainty 5.2% 

variation of track quality and geometrical cuts 0.5% 

choice of signal shape parameterization 2.0% 

error in number of B and B 1.8% 

choice of the 5m region widths 1.2% 

statistical error on efficiency (combined) 1.0% 

Total TU% 
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TABLE VI. Comparison of our results on "product branching ratios" (in %) with those of previous measurements. 
The error in this work's result is only the statistical one (see text). 



Channel 



Measured Br. Ratio 



This work 



B -> D U X 
B -> D+X 
B -> D*°X 
B -> D*+X 
B -> D*+X 



CLEO [15 1 



ARGUS Hi 



BbBd 
BbBd 
BbBd-Bd 
B b Bd>B d o 
BbB d *B d + 



2.51 ±0.06 
2.16 ±0.08 
0.620 ±0.031 
0.570 ±0.021 
0.639 ± 0.054 



2.33 ±0.12 ±0.14 
2.26 ±0.30 ±0.18 

0.56 ± 0.03 ± 0.05 



1.94 ±0.15 ±0.25 
2.09 ±0.27 ±0.40 

0.71 ±0.06 ±0.12 



TABLE VII. Inclusive B decay branching fractions to D and D*. 



B(B -> D°X) 
B(B -> D+X) 
B(B -> D*°X) 
B(B -» D*+X) 



0.636 ± 0.014 ± 0.019 ± 0.018 
0.235 ± 0.009 ± 0.009 ± 0.024 
0.247 ± 0.012 ± 0.018 ± 0.018 
0.239 ± 0.015 ± 0.014 ± 0.009 



TABLE VIII. Charm counting 



B{B -» D°X) + B(B w D + X) 
B(B -> D a X)S 
B(B -> A C X) & 
B(B H+X) I 
B(B - H°X) p 
2 x B(B -> Charmonia X) Q 



0.871 ± 0.035 
0.118 ±0.031 
0.039 ± 0.020 
0.008 ± 0.005 
0.012 ± 0.009 
0.054 ± 0.007 



2 We took the product of branching ratios measured by CLEO Q: B(B -> D S X)-B(D+ -» (jm + ) = (0.424±0.014± 
0.031)% and divided it by the world average B(D+ -» <t>ir + ) = (0.36 ± 0.09)% @. 

3 We have used the product branching ratio B(B A c X)B(A C -» pif ~tt+) = (0.181 ± 0.022 ± 0.024)% from Q 
and divided it by the world average B(A C -> pK~-K + ) = (4.4 ± 0.6)% Q. 

4 This branching fraction and the following one are derived from the product branching ratios reported in |p8[ . The 
B(B — > H c X) branching fractions in these papers were calculated under the assumptions that r(H c -^Xlv) = 
r(D — > Xlv). M.B. Voloshin |39| estimated that, because of Pauli interference of the strange quark, F(H C —* X£v) 
is in fact larger. His estimate (M.B. Voloshin, private communication) is T(H C — > Xlv) = (2.5 ± 1 ,0) r(D — > X Iv). 
Furthermore, here we assume r(S c — > "Elv) jY{E. c ~ > Xlv) — 0.8 ± 0.2 rather that 1.0 as in ref. |3q]. Altogether, 
the values of B(S + — » H _ 7r + 7r + ) and B(H° — > H _ 7r + ) are decreased by a factor of 2.0 ± 0.9 with respect to those in 
ref. |§. 

5 We took the inclusive, direct B decay branching fractions to ?/>/, %ci an d Xc-i from the CLEO result |^| and 
added 0.009 ± 0.003 (the upper limit on B(B — » r\ c X)) to take into account all other charmonium production from 
B decay. 
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FIGURES 
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FIG. 1. The m(K~TT + ) distribution for D° candidates (0.0 < xd < 0.50) (where xo = po/Pmax ) from "on 
resonance" data (p max = 4.950 GeV/c) and from "continuum" data (p max = 4.920 GeV/c). The "continuum" 
distribution is scaled by the luminosity and cross section factor 2.080. 
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FIG. 2. The 8m distribution for D* + D°ir + candidates (data). 
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FIG. 5. The "background-free" m(K~TT + ) distribution for D° mesons that are decay products of D* H 
0.0 < xd < 0.50. The two-Gaussian fit of the distribution is also shown. 
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FIG. 6. The "background-free" m(K~ tt + tt + ) distribution for D + mesons that are decay products of D* + — > D + tt" . 
0.0 < xd < 0.50. The two-Gaussian fit of the distribution is also shown. 
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FIG. 7. Representation of the functions used to fit the D signal and the various backgrounds for the momentum 
bin 0.050 < x < 0.075. Note the offset of the vertical scale. 
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FIG. 8. Representation of the functions used to fit the D° signal and the various backgrounds for the momentum 
bin 0.325 < x < 0.350. Note the offset of the vertical scale. 
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FIG. 9. The m(K 7r + ) distribution for D — ► K + tt after switching particle identity. Points with error bars are 
from data, the histogram is obtained from Monte Carlo using track tagging (MC-tag). Momentum bin 0.05 < x < 0.15. 
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FIG. 11. Switched if7r background. The ratio of the total area of the switched mass distributions (satisfying our 
particle identification selections) to the signal area, vs scaled momentum x = p/p ma x- 
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FIG. 12. The m(K~-K + ) distribution for B — ► D° X candidates in ten x bins from to 0.25. The line is the result 
of the fit. 
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FIG. 13. The m(K 7r + ) distribution for B — > D° X candidates in ten a; bins from 0.25 to 0.50. The line is the 
result of the fit. 
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FIG. 14. D° finding efficiency as a function of scaled momentum x. The histogram is the result of the smoothin 
fit, binned. 
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FIG. 15. The background subtracted and efficiency corrected yield of D° mesons from B 
function of scaled momentum (p m ax = 4.950 GeV/c). 
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FIG. 16. The m{K~ n + n + ) distribution for D + candidates from data at the T(45) resonance and from below BB 
threshold (scaled), 0.0 < xn < 0.5. The line is the result of the fit. 
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FIG. 17. The m(K~TY + Tv + ) distribution for B -> D + X candidates in ten x bins from 0.0 to 0.50. 
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FIG. 18. D + finding efficiency as a function of scaled momentum x. The histogram is the result of the smoothing 
fit, binned. 
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function of scaled momentum (p m ax = 4.950 GeV/c). 
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FIG. 20. D*° finding efficiency as a function of scaled momentum x. The histogram is the result of the smoothing 
fit, binned. 
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FIG. 21. The m(K~ 7r + ) distribution for D° candidates from B —* D*° X in the momentum interval 
(0.0 < x D * < 0.50). 
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FIG. 22. The m(K~n + ) distribution for D*° in ten x bins from to 0.50. The line is the result of the fit. 
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FIG. 23. The background subtracted and efficiency corrected yield of D*° mesons from B — ► D*°X — > (D°7r°) +X 
decay, as a function of scaled momentum (p max = 4.950 GeV/c). 
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D°7r + finding efficiency as a function of scaled momentum x. 
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FIG. 25. The m{K--n+) distribution for D° candidates from B -> D*+ X -> (D%+) X in 7 a; bins from 0.15 to 
0.50. The 8m sidebads have been subtracted. The line is the result of the fit. 



34 



6000 



-i 1 1 r 



153D1B7-D2B 
T 1 1 1 



5000 



4000 - 



in 
o 
o 



3000 



2000 



1000 



_l I l_ 



_l I I I I L. 

0.2S 
* =P'P m „ 



0.50 



FIG. 26. The background subtracted and efficiency corrected yield of D* mesons from B ~* D*+X -> (D°ty+) X 
decay, as a function of scaled momentum (p ma x = 4.950 GeV/c). 
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FIG. 27. The m{K~-K + ) distribution for D* + candidates from B decay (0.15 < xd* < 0.50), after 8m sideband 
subtraction. The dashed line shows the background under the signal. 
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FIG. 28. D*+ -> D+tt° finding efficiency as a function of the scaled momentum x. The histogram is the result of 
the smoothing fit, binned. 
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FIG. 29. The background subtracted and efficiency corrected yield of D* mesons from B -f D*+X (D+tt°) X 
decay, as a function of scaled momentum (p ma x = 4.950 GeV/c). 
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FIG. 30. The m(K tt + tt + ) distribution for D* + candidates from B decay (0.0 < xo* < 0.50), after &m sideband 
subtraction. 
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FIG. 31. The background subtracted and efficiency corrected yield of D* + mesons from B — > D* + X decay, as a 
function of scaled momentum (p max = 4.950 GeV/c). For < x < 0.15, only the D* + — > D + 7r° measurement is used. 
For 0.15 < x < 0.5, measurements from both D* + — > D°n + and D* + — * D + tt° are combined. 
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FIG. 32. The background subtracted and efficiency corrected yield of D* + as a function of cos for four x bins 
from 0.1 to 0.5. 
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FIG. 33. The background subtracted and efficiency corrected yield of D*° as a function of cos 6 for five x bins from 
0.1 to 0.5. 
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FIG. 34. a(D*) vs x = p/p max \Pmax = 4.950 GeV/c). 
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